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CALCULATED  PHOTON  SPECTRA  AND  RADIATION  DOSE 
FROM  50  MeV  ELECTRONS  STOPPING  IN  AIR 


I.  Introduction 

A Monte-Carlo  electron-photon  shower  calculation  has  been  performed 
with  the  computer  code  ETRAN  16  of  Berger  and  Seltzer to  determine 
the  bremsstrahlung  spectrum  and  dose  resulting  from  discharging  a 
monodirectlonal  beam  of  50  MeV  electrons  into  air.  The  purpose  of  the 
calculation  was  to  estimate  the  amount  of  acreage  and  approximate  shape 
of  an  accelerator  site  based  on  accepted  standards  of  radiation  safety. 
It  is  assumed  that  it  is  of  interest  to  Inject  the  beam  into  the  air 
outside  the  accelerator  building  and  because  it  is  not  feasible  to 
employ  conventional  radiation  shielding,  it  is  necessary  to  rely  on 
air  absorption  of  the  radiation  over  long  distances  to  reduce  the  dose 
to  harmless  levels.  Therefore,  the  minimum  area  occupied  by  the 
accelerator  site  is  an  Isodose  line,  corresponding  to  distances  from 
the  beam  which  meet  accepted  standards  for  safe  radiation  levels  given 
some  typical  electron  current  and  operating  conditions.  In  this  report 
we  have  used  the  non-radiation  worker  allowable  dose  level  of  0.5  Rem 
per  year. 

We  performed  the  electron-photon  shower  calculation  using  the  one 

dimensional  computer  code  ETRAN  16.  Because  this  code  Is  not  entirely 

suitable  for  a three  dimensional  system,  some  assumptions  must  be  made 

to  get  useful  results.  Namely,  we  have  assumed  that  since  we  are 

Interested  In  the  dose  at  distances  far  from  the  primary  photon  source, 

we  can  treat  the  photon  distribution  generated  after  one  electron  range 
Note;  Maniucript  submitted  July  30,  1976. 


as  a point  source  distribution  and  assume  simple  exponential  attenuation 
for  distances  further  away.  Because  no  previous  calculations  had  been 
performed  and  the  code  was  available  at  NRL,  we  felt  it  would  be  worth- 
while to  present  these  results,  which  can  serve  as  a good  estimate  of 
the  dose  at  0°  and  upper  limit  for  other  angles.  Our  use  of  this 
particular  version  is  described  in  Section  II  and  the  computational  re- 
sults from  ETRAN  are  presented  in  Section  III.  The  approximations  required 
to  apply  the  ETRAN  results  to  our  particular  geometry  and  the  conversion 
of  the  photon  spectra  to  dose  are  presented  in  Section  IV.  A dis- 
cussion of  results  is  presented  in  Section  V.  The  conclusions  and 
recommendations  are  given  in  Section  VI  and  Section  VII,  respectively. 

II.  Description  of  ETRAN  Computer  Code  Calculation 

The  computer  code  ETRAN  16  of  Berger  and  Seltzer was  used  to 

calculate  the  distribution  of  electrons  and  photons  emerging  from  a 

semi-infinite  slab  of  air  when  a monodlrectlonal  beam  of  50  MeV 

electrons  is  incident  normal  to  the  slab.  In  Fig.  1 the  electron  beam 

is  shown  incident  to  the  plane  at  Z « 0.  The  photons  are  shown  emerging 

from  the  solid  angle  defined  by  0 and  9 + A9  from  the  plane  at  Z ■ Z^. 

The  thickness  of  the  slab  Is  Z and  Infinite  In  the  other  two  dimensions. 

o 

In  using  ETRAN  various  options  in  the  codes  were  selected  to 
optimize  conditions  of  Interest  to  us.  A description  of  the  computer 
code  Including  these  various  options  have  been  given  In  detail  by  the 
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(1  2) 

authors  in  other  publications.  ’ Only  those  directly  reUvani 
to  our  calculation  will  be  mentioned  here. 

We  are  primarily  Interested  in  the  radiation  dose  produced  by  the 
photon  at  distances  much  greater  than  the  electron  range  in  air.  The 
dominant  contribution  to  the  radiation  dose  in  a phantom  at  these 
larger  distances  will  be  the  bremsstrahlung  radiation  emitted  as  the 
electrons  are  slowing  down  and  0.511  MeV  radiation  from  the  annihilation 
of  positron-electron  pairs.  Primary  electrons  do  not  directly  con- 
tribute to  the  radiation  dose  because  the  points  of  Interest  are  beyond 
their  range.  For  the  same  reason  knock-on  electrons  and  photon-produced 
electrons  do  not  directly  contribute.  In  ETRAN  the  Monte  Carlo  history 
of  each  primary  electron  was  determined  by  random  sampling  from  many 
interactions  taking  place  within  a short  path  traversed  by  the  electron 
instead  of  sampling  each  individual  interaction.  The  energy  loss  of 
the  electron  is  determined  by  using  a model  that  assumes  that  the 
electrons  are  slowing  down  continuously  rather  than  in  discrete  steps 
and  the  electron  angular  deflection  is  determined  by  sampling  from  the 
Goudsmlt-Saunderson  multiple-scattering  distribution.  The  histories  of 
the  bremsstrahlung  photons  produced  by  the  slowed-down  electrons  were 

generated  by  sampling  individual  interactions  according  to  known  formulas 

(3) 

of  continuous  bremsstrahlung.  The  photo-electrons,  pair  produced  particles. 


compcon  scattered  electrons,  and  knock-on  electrons  are  not  followed 
after  they  are  created.  But  the  bremsstrahlung  radiation  and 
annihilation  radiation  resulting  from  their  assumed  straight  line 
trajectories  are  included. 

Primary  electron  histories  are  terminated  when  the  electron  kinetic 
energy  falls  below  5.0  MeV  and  the  photon  histories  when  the  energy 
falls  below  0.200  MeV.  Consequently,  no  Auger  electron  of  K x-rays 
were  included  in  the  shower  calculation,  since  their  energies  were  below 
our  chosen  threshold  values.  It  is  estimated  that  the  exclusion  of 
particles  and  photons  or  x-rays  below  the  above  thresholds  and  the 
neglect  of  tracking  the  photon  produced  particles  would  affect 
the  dosage  by  much  less  than  50%. 

III.  Results  of  ETRAN  Calculation 

The  results  calculated  using  ETRAN  are  in  the  form  of  photon 
spectra  differentiated  in  photon  energy  and  angle.  The  emergent  angle 
of  the  photon  is  measured  with  respect  to  the  normal  to  the  plane 
surface  (see  Fig.  1).  Furthermore,  the  photon  spectrum  Includes  all 
photons  emerging  at  the  same  angle  summed  across  the  plane,  since 
ETRAN  16  does  not  keep  track  of  the  spatial  position  along  the  plane. 

In  Fig.  2 the  calculated  number  of  photons  per  MeV-sr  are  plotted 
versus  photon  energy  in  the  form  of  a histogram  for  several  different 
forward  angles  at  a distance  of  1 r^,  where  r^  is  the  mean  range  in 
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air  at  STP  for  50  MeV  electrons  (r^  150  meters).  At  each  given 

angle  above  some  threshold  energy,  the  number  of  photons  per  MeV-sr 
decreases  monotonically  as  the  photon  energy  Increases.  This  decrease 
becomes  more  rapid  as  the  angle  Increases.  These  histograms  also  show 
that  the  spectra  fall  off  very  sharply  as  the  angle  increases. 

In  Fig.  3 the  number  of  photons  per  sr  summed  over  energies  is 
plotted  versus  distance,  ^n  a given  angular  bin  the  points  fall  off 
in  distance  monotonically  and  can  be  fit  by  a straight  line  on  a semi- 
log plot.  This  is  reminiscent  of  simple  exponential  atenuation  of 
beam  of  photons  traversing  a thin  target. 

IV.  Calculation  of  Radiation  Dose  from  the  Photon  Spectra 

The  radiation  dose  absorbed  by  the  passage  of  the  gamma  rays  or 
photons  through  a small  phantom  sample  of  water  can  be  determined  by 
assuming  simple  exponential  absorption  using  known  mass-energy  absorption 

(4) 

coefficients.  The  dose  in  units  of  roentgen-equivalent  man  (Rem) 
per  electron  can  then  be  written  as 

D - E d(E  ,^,r)  AE 

E ' ’ 

Y 


where 


and  K ■ 1.60  X 10“®  gm/MeV. 

- energy  of  photon  bin  in  MeV, 
d(E^,!ji,r)  ■ number  of  photona/MeV-sr  at  the  point  p, 

AE^  ■ width  of  photon  energy  bin, 

u (E  ) ■ mass-energy  absorption  coefficient  in  water, 
w Y 

p - 1 gm/cm-*, 

s ■■  30  cm, 

r ■ distance  between  the  point  on  the  plane  surface 
P(r,$)  and  the  source, 

($  ~ angle  between  the  normal  to  the  surface  and  a line 
drawn  from  the  source  distribution  to  the  point  on 
the  surface  P(r,ij») 

The  biological  effectiveness  of  the  gamma  radiation  is  taken  as 
one  so  that  1 Rem  <■  1 Rad,  and  1 Rad  • 100  ergs/gram.  In  order  to 
determine  the  dose  at  some  point  in  space  it  is  necessary  to  know  the 
incident  photon  spectrum  N(E^,(|),r)  at  the  point  in  question.  The 
difficulty,  however,  is  that  ETRAN  16  calculates  N(E^,9,Z)  which 
includes  all  photons  emerging  at  angle  0 anywhere  along  the  plane 
(see  Fig.  1).  Clearly,  this  is  too  many  photons,  since  they  all  do 
not  contribute  to  the  dose  at  a given  point  on  Che  surface.  On  the 
other  hand  there  are  photons  that  are  emergent  at  other  angles  that 
should  be  Included  which  we  have  not  taken  into  account.  Most  of 


chese  phocons  do  not  originate  from  the  general  area  where  the  beam 
Is  actually  slowing  down  and  therefore  we  expect  Chese  contributions 
are  small  compared  to  chose  phocons  chat  cake  an  approximate  straight 
line  route. 

To  circumvent  this  difficulty  we  have  used  the  code  to  calculate 

the  photon  spectrum  emergent  at  an  angle  9 integrated  over  the  plane 

at  Z ■ Ir^,  a distance  at  which  all  primary  electrons  have  stopped. 

Since  we  are  Interested  in  the  dose  when  r >>  r^,  we  can  regard  the 

dimensions  of  the  region  150  meters)  in  which  the  electron  beam  is 

slowing  down  is  small  compared  to  r as  is  shown  in  Fig.  4.  We  now  use 

the  calculated  photon  distribution  as  a source  distribution,  which 

propagates  outwards  as  a point  source,  and  whose  intensity  decreases 

exponentially  with  distance  using  standard  air  mass-attenuation  co- 
(4) 

efficients.  Therefore,  the  photon  spectrum  at  P(r,i^)  can  be  written 
as 

N(E^,(|i,r)  a N(E^,9,r^)  exp  [-u^(E^)r]  . 

The  function  d(E^,(>,r)  becomes 

[-u  (E  )sl  [-U  (E  )rl 
d(E^.^.r)  - KE  N(E^.9,r^)(l-exp‘-  " jexp^  " " ■' 


where  N(E^,9,r^)  is  the  photon  spectrum  emergent  at  angle  9 and 
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integrated  over  the  plane  at  Z - Ir^  and  p^(E^)  is  the  gamma  ray 
mass-attenuation  coefficient  in  air.  We  expect  the  exponential 
attenuation  prediction  to  work  best  at  the  more  forward  angles  because 
here  the  gamma  rays  are  most  Intense  and  those  that  contribute  at 
these  angles  undergo  fewer  scatterings  than  those  at  the  wider  angles. 
The  dose  is  plotted  in  Fig.  5 for  angular  bins  of  0°-0.5°,  5°-10°  and 
25°-30°. 

The  horizontal  line  in  Fig.  5 represents  the  radiation  level  of 

0.5  Rem  per  year,  which  is  the  allowable  dose  a non-radlatlon  worker 

can  accumulate  over  a period  of  one  year.  This  line  is  calculated 
18 

assuming  that  10  electrons  per  day  are  injected  into  the  air  for  a 
period  of  one  year. 

V.  Discussion  of  Results 

The  intersection  of  the  horizontal  line  and  each  curve  gives  the 
distance  at  which  one  would  receive  the  allowable  dose.  A plot  of  this 
distance  as  a function  of  angle  is  shown  as  the  full  line  in  Fig.  6. 

The  Isodose  line  is  strongly  forward  peaked  at  0°  at  3300  meters  and 
falls  off  to  1000  meters  at  25*^.  Using  simply  exponential  attenuation 
excludes  y rays  scattered  into  the  angular  bin  or  so  called  "build- 
up". We  do  not  expect  this  to  be  too  important  between  0°  and  25°, 
since  the  primary  radiation  is  too  intense  in  these  directions.  How- 
ever, it  could  be  significant  at  wider  angles  where  the  primary 
radiation  is  much  weaker.  To  get  an  estimate  of  the  build-up  at 
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forward  angles,  we  carried  Che  MonCe-Carlo  calculation  out  Co  20 
ranges  and  compared  Che  results  Co  using  simple  exponential  attenuation. 
These  results  are  shown  In  Fig.  6 as  the  dashed  curve.  At  5°-10° 
the  more  complete  calculation  produced  a lower  dose  by  a factor  of  two 
and  at  25°-30°  a higher  dose  by  a factor  of  four  at  a distance  of  20 
ranges.  This  reduces  the  safe  distance  at  3°  by  57,  and  Increases  It  at 
25°  by  202.  We  Interpret  the  difference  at  25°  as  a breakdown  In  our 
model  assumptions  as  we  depart  from  the  forward  angles.  That  Is  Che 
"scattered  In"  contributions  are  becoming  greater  than  that  "scattered 
out";  Che  photons  In  general  are  undergoing  many  more  scatterings  so 
that  we  no  longer  have  approximate  straight  line  trajectories  between 
Che  source  and  the  point  In  question  Invalidating  the  point  source 
approximation.  The  dashed  curve  cannot  be  considered  completely  correct 
because  It  Includes  photons  Integrated  over  the  entire  seml-lnfinlte 
plane. 

At  forward  angles,  9 * 0°  to  25*^,  we  expect  our  calculated  dose  to 
be  fairly  realistic  and  considering  that  it  may  be  necessary  to  direct 
the  electron  beam  over  a wide  angular  range,  the  maximum  calculated 
safe  distance  of  3300  meters  is  certainly  a useful  upper  limit  for  the 
forward  directions.  The  forward  angles  may  be  used  as  upper  limits 
for  Che  backward  angles,  although  overestimating  the  acreage,  but  still 
safe  upper  limits. 
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VI.  Conclusions 


We  conclude  Chat  a distance  of  3300  meters  at  0*^  la  Indeed  a safe 
distance.  Since  the  safe  distance  would  be  expected  to  decrease 
rapidly  as  Che  angle  Increases,  as  the  results  show,  a distance  of  3300 
meters  can  be  used  as  a conservative  upper  limit  for  the  backward  angles. 
More  realistically  we  have  chosen  a distance  of  1000  meters  at  angles 
of  30°  and  greater. 

VII.  Recommendat Ions 

We  recommend  that  In  order  to  obtain  more  accurate  radiation  doses 
and  appropriate  safe  distances  at  angles  of  25°  and  greater  Chat  a 
complete  three  dimensional  Monte  Carlo  electron-photon  shower  calculation 
be  performed. 

For  radiation  doses  at  reflected  angles  (e.g.  0 > 90°)  and  for 
higher  bombarding  energies  (e.g.  > 100  MeV) , the  Inelastic  electron 

scattering  mechanisms  and  the  effect  of  the  finite  size  of  the  nucleus 
on  Che  elastic  electron  scattering  cross  section  becomes  progressively 
more  Important  as  9 and  Increase.  None  of  the  existing  electron 
transport  codes  Include  these  effects. 

At  NRL  we  are  In  the  process  of  using  an  electron  transport  code 
more  suitable  for  three  dimensional  cyllndrlcally  S3nnmetric  systems, 
and  also  considering  modifying  the  code  to  Include  Che  Inelastic 
scattering  and  the  effect  of  the  finite  size  of  the  nucleus. 
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Fig.  2 — A histogram  plot  of  the  number  of  photons/MeV-sr 
per  incident  electron  as  a function  of  photon  energy  for  var- 
ious angular  bins  integrated  over  a plane  at  Z =>  Ir,, 
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Fig.  4 — Target  geometry  for  the  small  source  approximation 
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Fig.  5 — A plot  of  the  calculated  radiation  dose  in  Rem/e  for 
different  angular  bins  as  a function  of  distance 
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Fig.  6 — A plot  of  the  “radiation-safe  distance”  as  a function  of  angle 
measured  from  the  incident  electron  beam  direction 


